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TECHNICAL MEMORANDUM X- 64646 


PASSIVE STABILITY OF A SPINNING SKYLAB 
INTRODUCTION 


In 1970 NASA’s Marshall Space Flight Center initiated a study to deter- 
mine the feasibility of spinning the Skylab (the first U.S. manned orbiting space 
station) . The purpose of the spin would be to provide artificial gravity so that 
its effects might be assessed. One facet of this study is presented in this paper 

In spinning the Skylab, it would be necessary to point the solar panels 
toward the sun. This would require the vehicle to spin about a principal axis 
of intermediate moment of inertia, which cannot be done stably. That axis can 
be made the axis of maximum moment of inertia by deploying masses on the 
tips of extendable booms in a direction perpendicular to the spin axis. A sim- 
plified model of the modified spinning Skylab is described, and analyses of the 
stability of motion and the rotational dynamics are presented. 


MODEL DEFINITION AND EQUATIONS OF MOTION 


The geometrically complex Skylab (j](Fig. l) is simplified to make it 
analytically tractable. The vehicle is modeled (Fig. 2) as a single rigid core 
body with principal moments of inertia I* (throughout this paper, subscript i 

ranges from 1 through 3 and refers to the 1, 2, 3 body-fixed axes) with 
Ii <I 3 <I 2 . Attached to the core body are two flexible massless booms, each of 
length T (measured along the 2-axis in steady state) and each having a tip mass 
m. In the steady state, the principal moments of inertia of the entire vehicle 
coincide with the body-fixed axes and are designated I .with IpcI^I^ The vehi- 
cle is assumed to have a steady-state spin velocity fi about the 3-axis, and the 
angular velocity vector is written in body-fixed coordinates as co = (w t , w 2 , fi 
+ w 3 ), where w. represents small perturbations about the steady state, i. e. , 

t k 

|wj«l. Displacements u, (k = 1,2) of the tip masses from the spinning 

steady state are assumed to be small. Nonrotating boom stiffness is charac- 
terized by stiffness coefficients k.. Structural damping is assumed proportional 
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Figure 1. Skylab. 



Figure 2. 


Simplified model. 





to elastic deformation velocities and is denoted by coefficients d.. Rotational 

1 

dynamics of the vehicle may be represented by a set of nine differential equa- 

k 

tions written in variables ir , w^. The set may be reduced to six by making the 

substitution u. = u. 1 - u. 2 , where u. represents the skew symmetric mode of the 
1111 

elastic deformations and hence causes angular motion about the vehicle's 
steady-state attitude. In this investigation stability of rotational motion is of 
interest, so only the skew symmetric mode will be considered. The corre- 
sponding linearized equations of motion are 

wj - K t w 2 + y jOf/iy + g 3 ) = Ti/Ijft 

w<l - KjWj = T 2 /I 2 fl 

(l/fi)(w| + w 2 ) + n'j + + (cr 3 2 + l)p 3 =0 (1) 


and 

w 3 - Vs 0 (mV - ^g 2 ) = t 3 /i 3 J2 

(-l/fi)w^ + n\' + A + c T?Hi - (1/r 2 )g 2 ' = 0 

( -2/fi) w 3 + (l/r 2 )gj + jt' 2 ' + A 2 p^ + (<x 2 - l)p 2 = 0 , (2) 

where 

I t = if + 2mT 2 
12 = 12 * 

I 3 = if + 2mT 2 

cr 2 = k./ mQ 2 

i i 

y. = 2mr 2 /I. 

ii 

A. = d./mf2 
i i 

K t = (I 2 - Ij)/!! 
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K 2 = (I 3 - Ii)/I 2 


ii. = u./2 r 
1 1 

and T represents the applied torques. For physical reasons, 
i 

|K k l< i (k = 1,2), 

0 < y < 1 (m = 1, 3), 
m 

d. > 0, 


and 


k. > 0. 

1 

Primes represent derivatives with respect to r, where 
r = fit, 

and t denotes real time. Equations (l) describe wobble motion (resulting from 
angular motions about the 1- and 2-axes and linear motion in the 3-direction) . 
Equations (2) describe in-plane motion (a combination of angular motion about 
the 3-axis and linear motions in the 1- and 2-directions). 


STABILITY ANALYSIS 


The motion of the vehicle can be described by a nutation about the axis 
of angular momentum which, in the absence of external torques, is inertially 
fixed. This motion is called passively stable if the nutation damps out and the 
vehicle rotates only about the axis of angular momentum (the effect of applied 
torques T^ is to change the attitude of this axis). The following analysis 

determines under what conditions the spinning vehicle can be passively stabi- 
lized, assuming I 1 <I 2 <I 3 . 
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The stability of wobble motion is investigated by obtaining the charac- 
teristic equation from equations ( l) : 

(1 - Tl )* 4 + A 3 \ 3 - (K t K 2 + Yl K 2 + y x - a l - l)X 2 - AjK^X 

-[(a 3 2 + lJjiqKj + y t K* ] « 0, (3) 

where 

X = s/S2 = 7} + iy 

and s is the Laplace operator. Regions of stability may be determined by 
applying the D-composition technique f 2] . The stability boundary associated 
with the origin of the complex X-plane is found by setting the coefficient of X° 
equal to zero in equation (3), yielding 

K, = -[ (o-3 2 + lVyjKiKj. (4) 

The boundary associated with infinity on the X-plane is found by setting the 
coefficient of X 4 equal to zero: 

7l = 1. (5) 

The stability boundary associated with the imaginary axis of the X-plane is 
found by setting X equal to iy in equation (3) and writing the real and imaginary 
parts of the result as the two equations 

[y 2 - (cr 3 + 1 )] K t K 2 + Ti(y 2 - l)K 2 = -v 2 [( 1 - y Jy 2 + ( y t - a 3 2 - l)] 

and 

-AjyfoKg = -y 2 ). (6) 

Equations (6) may be solved for the relationship 

K 2 = -K,K 2 . (7) 

Finally, a stability boundary may exist that is associated with the singular case 
where equations (6) are not independent; i. e. , when the Jacobian J becomes 
identically equal to zero [ 2] . The Jacobian of equations ( 6) is 

J = A 3 y 1 y(y 2 - l), (8) 
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which is zero when A 3 , y it or p equals zero or when p 2 = 1. If p 2 = 1 is 
substituted into equations (6), the J = 0 stability boundary is found as 

K t K 2 = -1. (9) 

These stability results corroborate those reported by Barbera[3] . 

The wobble motion stability requirements may be portrayed on a-K 1 K 2 , 

K 2 parameter plane (Fig. 3). The stable region may be found by considering 
the sign of J [2] . Double cross-hatching is used on the boundary associated 
with the imaginary axis (including the J = 0 boundary) „ The hatching lies on 
the side of the boundary that is toward the stable region; in the X-plane, it lies 
on the left side of the boundary for increasing p. If J>0, the hatching in the 
parameter plane will also lie on the left side of the boundary as p increases 
(and on the right side if J<0). Single cross-hatching is used similarly on the 
boundary associated with the origin of the \-plane. Continuity of cross-hatching 



1. CROSS-HATCHED REGION IS STABILITY REGION 
LOST BECAUSE OF FLEXIBILITY 

2. ENCIRCLED NUMBERS REFER TO NUMBER OF STABLE 
ROOTS OF CHARACTERISTIC EQUATION 

3. ADDITIONAL CONSTRAINTS: 0<y ] <l,A 3 >0. 

Figure 3. Stability region for wobble motion. 
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must exist at contour intersections in the parameter plane corresponding to the 
origin of the A. -plane and to other intersections where the contours have equal 
values of v . 

The effect of boom flexibility on stability may be assessed by consider- 
ing the booms to be rigid; i. e. , setting l/o\ 2 equal- to zero. The boundary of 

equation (4) becomes = 0. The stable region is thus decreased by boom 
flexibility, as shown by the shaded region of Figure 3. Also, the effect of boom 
length and tip mass may be seen by considering their effect on K lt 

Similarly, stability of the in-plane motion may be determined by examin- 
ing the characteristic equation obtained from equations (2). 

A{ ( 1 - y 3 )x 4 + [Aj + A 2 (l - y 3 )]\ 3 

+ [of + o- 2 2 (l - y 3 ) + 3(1 - y 3 ) + A^U 2 

+ [A t (of - 1 + 4y 3 ) + A 2 of]X + of (of - 1 + 4y 3 )} = 0. (10) 

If it is assumed that the stiffness coefficient k 2 is a times as great as k 1 and 
k 3 , where a is a positive constant, and that kj *= k 3 , the application of 
D-decomposition to equation (10) yields the boundary associated with the origin 
of the X- plane: 

of = 0 and of = 1 - 4y 3 . (ll) 

The boundary associated with infinity on the X-plane is y 3 = 1. The boundary 
associated with the imaginary axis of the X-plane yields the same two bound- 
aries plus an additional constraint, Ai = 0. No additional boundaries arise 
from the J = 0 condition. The resulting region of stability is shown on a y 3 , 
erf parameter plane (Fig. 4), which shows the effects of the boom characteris- 
tics on the stability of the in-plane motion. 


WOBBLE MOTION DYNAMICS 


Once stability has been assured, the transient dynamics of the wobble 
motion are of interest. The character of this motion may be investigated by 
using the parameter plane technique [2]. Characteristic equation (3) is used, 
selecting inertia ratios and K 2 as two adjustable parameters of interest to 
the designer; i. e. , the effect of various values of KjKj and K 2 on the wobble 
motion dynamics will be determined. This is done by finding the relationship 
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NOTESs 

1. SEE NOTE 2, FIG. 3 

2. ADDITIONAL CONSTRAINT: A 1 > 0 

Figure 4. Stability region for in-plane motion. 


between these two parameters and the roots of equation (3). The relationship 
may be shown by mapping a convenient contour from the A- plane onto the KjK 2 , 
K 2 parameter pla,ne. The contour selected is the line 

A = + i„ n (l - £ 2 ) /z , (12) 


where v 

n 


v(l - £ 2 ) 



Any root of equation (3) lying on this contour will 


have a specified damping ratio f . This line is mapped onto the parameter 
plane by substituting equation (12) into equation (3), which may be written in 
the form 



(a K.K 2 + b K, + 

' p 1 i p i 




where values for coefficients a , b , c are found in Table 1. 

P P P 


(13) 
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TABLE 1. COEFFICIENTS OF WOBBLE CHARACTERISTIC EQUATION 


p 

a 

P 

b 

P 

c 

P 

0 

-(°3 + !) 

-Yi 

0 

1 

_a 3 

0 

0 

2 

-1 

-Yl 

1 + cr 3 2 - Ti 

3 

0 

0 

^3 

4 

0 

0 

1 - Yl 


If \ P is written as the sum of its real and imaginary parts, 

A P = X + iY , (14) 

P P 


values of X and Y may be found from the recurrence formulas, 
P P 

X ^ + 2£v X + v 2 X =0, 
p+1 n p n p-1 


Y _ + Y + v 2 Y =0, 
p+1 n p n p-1 


(15) 


where it is seen from equations (12) and (14) that X 0 = 1, Y 0 = 0, X t = , 

V n 

and Y t « (l - f 2 ) 2 . Substituting equation (14) into equation (13) and 

separating the resulting real and imaginary parts, one obtains two simultaneous 
algebraic equations that may be solved for parameters K t K 2 and K 2 : 


KiKj = (B 1 C 2 - B 2 C 1 )/J, K 2 = (A 2 Cj - A 1 C 2 )/J, 

j = a 1 b 2 — a 2 b 1? 


where 


A, = 


4 

7 , a x , 

a p p 


B > = 2 W 

p=0 F 


- 2 


c X , 

p p 

p=0 


(16) 
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and 


A 


2 


E 


p=0 


a Y , 
P P 


4 

- E 

p=0 


b Y , 
P P 


c 2 


4 


E 

P=0 


c Y . 
P P 


(17) 


Equations (16) define the parameter plane contour corresponding to equation 

( 12) in the \-plane for a chosen value of f and coefficients a , b , c . Each 
v P P P 

contour is plotted as a function of the nondimensional natural frequency v 


n 


The real roots of the system may also be mapped from the X-plane onto 
the parameter plane. Each selected value r\ of a real root is substituted for X 
in characteristic equation (T3). The resulting equation will define a contour on 
the parameter plane representing the selected location rj of the real root. 

The response of the system in terms of variables Wj, w 2 , may be 
found by casting equations (l) in the form 


x f = Fx , 

where the vector x is defined as 

T 

x = [w t , w 2 , n 3 , /ij] , 
and the matrix F is 

0 (Yi + K i) 7i<?3 TlA 3 

K 2 (l - 7i) 0 0 0 

0 0 0 (1 - 7i) 

0 -(Ki + 1) ( Tl - (r 3 2 - 1) -A s 
The response x(r) may then be found from 
x(t) = 0(t)x(o) , 


F = 


(i - n) 


(18) 


(19) 


( 20 ) 


( 21 ) 


where x(o) is the initial condition vector and 4>(t) is the state transition matrix 
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( 22 ) 


0(t) = e FT =£ -1 {[\E - F]" 1 }. 

Matrix E is a unit matrix of the same dimensions as the F matrix, and £ -1 
represents the inverse Laplace operation. In evaluating <P(t), the roots of 
the characteristic equation (3) are the roots of the determinant | XE - F | and 
may be taken directly from the parameter plane. For example, if initial con- 
ditions exist only on w 3 and w 2 ; i. e. , 

x(0) - (x 10 , x 20 , 0, 0) T , (23) 

the w 3 (t ) response is found from 

w t (r) = £ _1 [$n(X)x 10 + $ 21 (X)x 20 ], (24) 

where 

* u (\) = X[X 2 + A 3 X/(1 - y t ) + (1 + ai - y t )/(l - Ti)]/|XE - F|, 
$ 21 (X) = K 2 $ n (\)/\ . (25) 

Example 

■' The model described by equations (l) and (2) may be used to crudely 

~ represent a spinning Skylab. Representative values for the coefficients in 
equation (3) are given in Table 2. 

TABLE 2. PHYSICAL CHARACTERISTICS OF SKYLAB 


If = 1. 01 x 10 6 kg m 2 
I* = 6. 90 x 10 6 kg m 2 

I* = 6. 85 x 10 6 kg m 2 


T = 23. 3 m 
m = 227 kg 

kj = k 3 = 146 N/m 


k 2 = 7.4 x 10 4 N/m 
= 0. 04 (kjn) ^ 

= 0. 6 s -1 


If the values are substituted into equations (16), ^-contours may be 
plotted on the K 2 parameter plane as functions of v (Fig. 5). Stability 

boundaries defined by equations (4), (7), and (9) are also plotted in Figure 5. 
Real root boundaries are found to be outside the stable region and hence are of 
no interest in the Skylab problem. 
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Figure 5. K 2 versus K t K 2 parameter plane plot. 
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A typical design point for Skylab is shown on the parameter plane. It 
corresponds to the physical characteristics shown in Table 2. The correspond- 
ing roots of the characteristic equation are seen to have damping ratios (£) and 
nondimens ional natural frequencies of 


£ = 0. 002, v = 0. 2854 
n 

and 

t= 0.01863, v n = 1.7612. (26) 

If the values of Table 2 and equations (26) are substituted into equation 
(21), the response for the Skylab model is found. For example, from equation 
(24), 

w t (r) = 1. 0978 e t> °° 5672T cos(0. 2854r + 0. 3274)w 1 (0) 

+ 0. 09359 e _0 ‘° 328lT cos(l.76lT + 4. 1270)w 2 (0) , (27) 

where x(o) is defined by equation (23). For w t (o) = w 2 (0) = 0. 001, the 
response w^r) is plotted in Figure 6. Similarly w 2 (r) and ^ 3 (t) are also 
determined and shown in Figures 7 and 8. 

CONCLUSION 


It has been shown that it is possible to passively stabilize the motion of 
a simplified model of a spinning Skylab by deploying flexible booms, thus 
altering the moments of inertia. Analytical results indicate the required boom 
stiffness properties for given vehicle mass properties and spin rates to achieve 
passive stability. Further, the use of the simplified model leads to results 
amenable to physical interpretation. 

To gain confidence that these results will apply to the actual Skylab, an 
additional step is being implemented. A detailed digital simulation model of 
the spinning Skylab vehicle has been developed at the Marshall Space Flight 
Center, and results compare favorably with those of the simplified model. 
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